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A perturbation expansion is considered about the Ohta-Jasnow-Kawasaki theory of phase- 
ordering dynamics; the non-hnear terms neglected in the OJK calculation are reinstated and treated 
as a perturbation to the linearised equation. The first order correction term to the pair correlation 
function is calculated in the large-d limit and found to be of order 
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I. INTRODUCTION 



I When a system is quenched from a high-temperature, homogeneous phase, into a two-phase region, domains of the 
equilibrium phases form and evolve with time. If the order parameter is non-conserved, then the coarsening dynamics 
[ are modelled by the time-dependent Ginzburg-Landau (TDGL) equation 

I In one dimension this system is exactly soluble , but due to the non-linear nature of the TDGL equation, no 
' exact solutions are available for a general number of dimensions so we must rely on approximate theories. There 
^ , are several approximate theories which describe the coarsening dynamics of this system. However, they all rely on a 
fS| ' similar approach in which the order parameter, (j), is replaced by a smoothly varying auxiliary field, m, which has the 
[ same sign as (j), and whose zeros define the domain walls. The equation of motion (TDGL) is then recast in terms 
of this auxiliary field and various approximations are made in order to make the resulting equation soluble. In the 
Ohta-Jasnow-Kawasaki (OJK) |^ theory, this is achieved by replacing the non-linear terms by their spherical average, 
thus linearising the equation for the auxiliary field . 

Comparison of the OJK results with simulation data have shown that this theory gives a very good description 
^ of the system. However, the approximation is uncontrolled and it is unclear why this approach gives such accurate 

d predictions. Furthermore, Blundell et. al. |^] have pointed out that by plotting (1 — C^) against 1/(7^2, where 

C^ = (0(x + r,t)0(x,<)), (1) 

^ ((l-0^(x + r,t))(l-0^(x,t))) 
O : (l-02(x + r,t))(l-<^2(x^i)) ' ^ ) 

o ■ . . 

• • ■ thereby removing any adjustable parameters, the OJK theory does not follow the simulation data as closely as previous 
comparisons, using only C^, suggested. 

An obvious way to examine this approximation is to treat the neglected non-linear terms as a perturbation to the 
■ linearized equation of motion, and this is exactly the approach taken in this paper. 

This calculation has two main advantages: firstly, in OJK theory the initial conditions are conventionally taken to 
be Gaussian; since the evolution equation is linear, this will ensure that the distribution for the auxiliary field at later 
times is also Gaussian. This feature of OJK theory has been questioned by Yeung et. al. Their simulations, which 
explicitly calculate the auxiliary field distribution, give results which are not exactly Gaussian, particularly at small 
values of the auxiliary field. We note that recently, in a series of papers, Mazenko HJl^, and Mazcnko and Wickham 
PI , have presented an approximate theory which goes beyond the Gaussian distribution. 

Secondly, it has been proposed by both Bray and Humayun iQ, and Liu and Mazenko that the OJK approx- 
imation becomes exact as the number of dimensions approaches infinity. Evaluating the dimensional dependence of 
the first order correction term enables this hypothesis to be tested. 

The main result of this paper is that the first order correction to the correlation function is 0{l/cP), lending weight 
to the assertion that the OJK theory becomes exact in an infinite number of dimensions. The exact form of the 
leading order correction to the two-time correlation function, Ci(r,ti,t2), is calculated in the large-c? limit; this is 
found to obey both Porod's Law iQ, and the Tomita sum rule with the singular contribution in the correction 

term modifying the amplitude of the Porod tail, and in the limit ti ^ <2, it is found that the correction term is of 
exactly the same form as the zero order result. The complexity of the calculation, however, has so far prevented the 
evaluation of higher order terms. 
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In the following section we present an outline of the OJK calculation, which provides the starting point for the 
perturbation calculation, we then in section III proceed with a detailed description of the perturbation calculation. 
We conclude with a summary and discussion of the results. 



II. THE OJK THEORY 



Consider a system described by a non-conserved scalar order parameter. The evolution of this system following 
a rapid quench from a high temperature homogeneous phase to a regime where there are two equilibrium phases is 
described by the TDGL equation, 

where F[(j)\ is a Ginzburg-Landau free energy functional, given by 



J d^xQ|V</.|^ + y(0) 



(4) 



V{(j)) is a potential whose minima define the equilibrium values of the order parameter; the conventional choice in the 
case of a scalar order parameter is V{(f>) = (1 — 0^)^/4. Note that there is no noise term in this equation, hence the 
results are only valid for quenches to temperatures where the effect of thermal fluctuations is negligible. 

Following the work of Ohta, Jasnow and Kawasaki jsj, the scalar field, cj){v^t), is replaced by a smoothly varying 
auxiliary field, m(r, t), where (/)(r, t) and m(r, t) have the same sign. The zeros of the field m then define the positions 
of the domain walls, and the normal, n, to a domain wall in the direction of increasing (j) is given by 



Vm 



(5) 



Inserting this into the Allen-Cahn equation for the velocity of an interface |]T^, v = —K = — V.fi; K is the curvature 
of the interface, we obtain: 

1 / _2 m \ 

where v is the speed of the interface in the direction of n. We now seek to obtain an evolution equation for m(r, t) 
by linking the domain-wall velocity to the time-dependence of the auxiliary field. The total time derivative of the 
auxiliary field in a frame of reference moving at a velocity v is given by 

dm dm , , 

If this frame is moving with the interface velocity, then v — vn and Vm are parallel, and the total derivative of m 
vanishes, implying that 

dm(r,t) , 

' =-i^|Vm|. (8) 

We now substitute the expression for the interface velocity (from equation (H)) into equation to obtain, 

dm 2 d^m , , 

This is the time evolution equation for the auxiliary field obtained by OJK. 

To make analytic progress, we now linearise this partial differential equation by replacing Uirij by its spherical 
average Sij/d. In this approximation, the evolution of the field is governed by a simple diffusion equation with 
diffusion constant D — 1 — 1/d. To calculate the correlation function we need to express the original order parameter 
(f) in terms of the auxiliary field. In the thin wall limit this is given hy (j) — sgn(m), therefore 

C(r,ti,i2) = (sgn(m(x + r, ii))sgn(m(x, i2)))- (10) 
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It is convenient to choose Gaussian initial conditions for the auxihary field, with zero mean, and correlator 

(m(x + r,0)TO(x,0)) = A(5(r). (11) 

Since the evolution equation is linear, the auxiliary field will then have a Gaussian distribution at all times. Therefore 
to calculate the correlation function we simply need to evaluate the joint probability distribution for m(x+r, ti ) = m(l) 
and m(x, ^2) = m{2). This is given by 



P(m(l),m(2)) : 
X exp 

where 7 is the normalized correlator 
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2^(5(1)5(2)(1-7^))V2 

1 /m(l)2 m(2)2 



27- 



m(l)m(2) 



2(1-72) V 5(1) 3(2) '(5(l)5(2))i/2 
(m(l)m(2)) 



m(l)2)i/2(m(2)2)i/2 

d/4 

exp 



{ti + t2y 



^D{ti + t2) 



(12) 

(13) 
(14) 



5(1) = (m(l)2) and S{2) ^ {m{2f) §. Completion of the average over the field m yields the final OJK result for 
the correlation function, and is given by 



C(r,ti,i2) = - sin ^7. 



(15) 



We note that this result only has a trivial dependence on the dimension d, through the diffusion constant, D = 1 — 1 / d. 



III. PERTURBATION THEORY 

The starting point for this calculation is the OJK equation for the evolution of the auxiliary field to(x, t) (derived 
in section ||), 

dm 2 d'^m 

where n = Vm/| Vrnj. As noted in the previous section, this equation is highly non-linear and to make analytic 
progress OJK approximated the non- linear term, niUj, by its spherical average, Sij/d, reducing equation ( p6[ ) to the 
diffusion equation with diffusion constant D = 1 — 1/d. 

In this calculation, the terms dropped in this approximation are treated as a perturbation to the diffusion equation 
and the first order correction to the correlation function is calculated. Therefore the equation we wish to solve is 



dm 9 , / Sii \ d^m , , 



where A is a small perturbation parameter. 

The solution of equation ([l^ ) can be expressed as a power series in A, m(r, i) — tuq + Ami + 0(A2); substituting 
this back into equation ( p^ ) and equating orders of A gives two coupled partial differential equations for mg and mi, 
the solutions to which are given by 

mo(r, t)^ j d'^fi G(r - /x, t) nioifi, 0), (18) 
mi (r, t) = - y dr J d'^^ G(r - n,t- t) 



X 



dimo{fJ,,T)djmo{fi,T) Sij \ d'^mnifi^r) ^^^^ 



|Vmo(/x,T)|2 d J dfj.^dfj. 
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where ~ r,T ~ t) is the Green's function for the diffusion equation and is given by 



G(r,i) 



■ exp 



r 



(20) 



and dimo{fi,T) = ^"g^ff ■ 

In the thin wall limit the order parameter is related to the auxiliary field by the equation (f) — sgn(m). Hence, by 
using the integral representation of sgn(m) and expanding the expression for the correlation function in A, we find 
that 



where 



C(r,ii,t2)= (<^(x + r,ti)(/.(x,i2)) =Co(r,ii,f2) + ACi(r,ti,i2) + 0(A2), 

Co{r,ti,t2) = (sgn(mo(x + r,ti))sgn(TOo(x,t2))), 
Ci(r,ti,i2) = 2( (sgn(TOo(x + r, ii))(5(TOo(x,t2))"T-i(x, i2)) 



(24) 



(21) 



(22) 



+ (sgn(mo(x,t2))i5(TOo(x + r, ii))mi(x + r, ii)) ). (23) 

If we define Ci{r, ti,t2) by 

C'i(r, ti,t2) = 2(sgn(mo(x + r, ti)) (5(mo(x, 12)) toi(x, t2)), 

then the first order correction to the correlation function is given by 

Ci(r,ii,t2) = C'i(r,ti,i2) + Ci(-r,t2,ti). (25) 

Since the two terms on the RHS of the expression for Ci(r,ti,t2) (equation (p5|)) differ only in that ii ^2 and 
r — r, we will only deal with one of these averaged terms, Ci{r,ti,t2), evaluating the complete expression at the 
end of the calculation. 

Substituting for mi from equation ( |l^ into equation (p^, we obtain the following expression for Ci{r,ti,t2), 

C'i(r,ti,i2) = -2^ dr J d'^ fj. G{x - n,t2 - T){sgn[mo{x + r,ti)] 

'dimo{fi,T)d.jmo{n,T) SiA d'^mo{fi,T) 
X I r= — J I d[mo{x,t2))- 



(26) 



|Vmo(/x,T)p d J ' dfii,dfij 

We now impose the conventional Gaussian initial conditions, with a zero mean, and a correlator given by 

(m(x + r,0)TO(x,0)) = A(5(r), (27) 

— — from the 



The average on the RHS of equation ( pq ) can be evaluated if we extract the differential operator q^.q^. 
average by defining a new spatial variable u. Equation (Eq) can then be written as 



C'i(r,ti,t2) = -2 dr ( d^ii G(x-M,i2-T) ^ 
Jo J dvidvj 



where 



dtma{p,,T)djmo{fi,T) _ % 
|Vmo(At,r)|2 " d 



sgn(mo(x + r, h)) mo{v, r) (5(too (x, <2)) \ 



(28) 



(29) 



The ensemble average in equation ( p9| ) can now be completed by using the joint probability distribution connecting 
the variables mo(x + r,ti), too(x, i2), ^0(1^, r) and 



For simplicity, at this point we introduce a contracted notation: 
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'mo{i^,T) m(2), 
mo(x,t2) ^ m(3), 



mo(x + r,<i) ^ m(l), 
and define the vector m by 

m = (m'j^ , m'2 , TO3 , • • • , , m(l) , m(2) , m(3) ) . 
The joint probability distribution for the components of this vector is then given by, 

P(m) = 5-^ exp(--TOi Aij mj), 

(2^)^(detA-i)3 '2 



where A^^ = {rriim 



where 



■j). P(m) is evaluated in appendix ^ and is given by 



P(m) 



(27r)'^(detA-i)* 



- exp 



F(m) 



2detyl-i 



(30) 
(31) 

(32) 
(33) 

(34) 



F{m) = m'^ekim'i + 2|.m'r7i(l) + 2r].m'm{2) + 2C .m' m{3) + pm{l)'^ + qm{2f 
+rm{3f + 2sm{l)m{2) + 2tm{l)m{3) + 2um{2)m{3). 

and the expressions for the coefficients ^, t], ^ and p, q, . . . u are given in appendix 
Equation (E3) may therefore be written in the form 



(35) 



xP(m) 



dm(l) / dm{2) / rfm(3) 



) sgn[m(l)]m(2)<5[m(3)]. 



(36) 



The function C'i(r,ti,t2) defined by equation ( p8| ) is evaluated in three main steps. Firstly, the expression for /y 



(equation (|36|)) is evaluated. This result is then differentiated to obtain , and finally, this expression is 

substituted back into equation ( |2^ ) and the remaining integrals are completed. 

To evaluate equation (36), we must complete the integrals over rh. The integral over m(3) is trivial; on setting m(3) 
to zero in equation (|3^), we see that the integral over m(2) can be reduced to a Gaussian by completing the square 
in the exponent of the probability density function, —F{m.)/2detA~^. For simplicity we deal with the expression for 
F(m) (equation (^5|)) first. Completing the square using the substitution m'(2) = m(2) + (sm(l) + rj.m')/^ in m(2) 
gives 



F(m) = m'fc 0ki 



VkVl 



, , {Pq - S^) ,,x2 ^ (ST?- g|) , ,,„s2 

1; H to(1) — 2 .m m(l) + qm (2) . 



(37) 



This may be simplified further by substituting for ^ and r) in the factor sr) — g^, from equations (A27) and (A2S) 
respectively, to give 



ST] — q^ = {pq — s^)a — {us — qt)c = {pq ~ s^)h. 



where this equation is used to define the vector h. Using equations ( |A42| ) and ( A44 ), h may then be written as 



h = a 



( us — qt 
\pq — 



zy — a.c 

zXs - c2 



On substituting equation (BSh back into equation (p7|), P(rii) reduces to 



F{rh) = m'fc 9ki - 



VkVl 



{pq~ s^) 



i{lf - 2h.m' m(l) + qm'{2) 



(38) 



(39) 



(40) 



Substituting equations (34) and (M) back into equation (M) gives 
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-(d+3)/2 



(det.4-i)i/2 

X / dm(l) sgn[m(l)] exp 



2detA 



Ym^.{0ki )mi 



(pq - s2) 



X f dm'{2) (m'(2) 



2q det A 
(sm(l) + 17. m') 



^(m(l)2-2h.m'TO(l)) 





■ (jm'(2)2" 


j exp 


2detA-i 



(41) 



The m(2) integral is now in the form of a Gaussian, on completion of which, lij simplifies to 
-(27r)-(''+2)/2 



.3/2 



m' 



exp 



2detA 



—mk{Oki )m; 



X / o?to(1) sgn[m(l)](sm(l) + ?7.m') 
{pq - s^) 



X exp 



2q Act A 



^{m{\Y -2h.m'TO(l)) 



(42) 



We now manipulate this expression into a form in which the m' integral can be completed. Firstly we need to 
remove the |m'p term from the denominator; this can be achieved by rewriting the denominator as an integral over 
an exponential, 



1 

]^ = J dvexp[-v\in'\^] 



(43) 



The second step is to manipulate the m{l) dependence so that the exponent in the complete expression for lij can 
be written in the form m'^flm', where J7 is a d x d matrix to be determined later. 

Due to the presence of the sgn[m(l)] function, the positive and negative ranges of the integral must be treated 
separately; in each case we complete the square in the exponent using the substitution u = (m(l) — h.m')/h.m'. 
The anti-symmetric contribution to the m(l) integral in equation ( ^ ) may be evaluated directly, and the symmetric 
contribution reduces to a term proportional to an error function. This integral, therefore, reduces to 



dm{l) sgn[TO(l)](sm(l) + rj.m') exp 



(pg - 

2q detA- 



^(m(l)2-2h.m'm(l)) 



= a + jSkiTn'i^m'i / duexp 



{pq - s^) 



2gdet^- 



where 



a 



2qs det^-i 
{pq - s2) 



and (3ki = 2/1^(77/ + shi). 



(44) 



(45) 



Substituting equations (43) and ( p4| ) back into the expression for (equation ([4^)), we obtain 

1 



1 



(27r)(^'+2)/253/2 

X I - m^nij 



d m' exp 



2detA 



—m^{&ef )m 



q 



X a + Pki ^^k'^^i / dii exp 



J dwexp [— u|m'p]^ 
{pq-s^) 



^(h.m')nw -1) 



2qdetA- 

We are now in a position to define the general matrix 17 (u, v) referred to earlier; 



r2y({2,w) 



1 



detA- 



h,hj{u^ - 1) +2M 



(46) 



(47) 



By expanding the RHS of equation (^6|), we see that the exponent of each contribution is given by the general dx d 
matrix Q,{u,v), evaluated at different values of u and v. Equation ( |4^ ) can therefore be written as the sum of four 
standard integrals. 
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(48) 



where 



I? 



TV 

^2 



^3 



TV 

^4 



d(27r)('i+2)/2q3/2 



d'^m' exp 



-im'^17(l,0)m' 



27rdq3/2(detl7(l,0))i/2' 
a 



(27r)(d+2)/2g3/2 

a 



dv 



27rg3/2 7„ (detf^(l,5))i/2' 



^(27r)(d+2)/2g3/2 



du I d'^m' m'f.m'i exp 



-im''^17(l,{})m' 



-im'^17(M,0)m' 



2Trdqy^ Jo (detf}(M,0))i/2' 



(49) 



(50) 



(51) 



(27r)(''+2)/2g3/2 
Pkl 



m^mjm^.mi exp 



1 T 

-m' r2({t,-D)m' 



27rg3/2 7p 



du / du 



(detf7(u,{}))i/2 



(52) 



Details of the calculation of the inverse and determinant of Vl{u, v) are supplied in appendix They are given by: 



u kj k 



in, J 



(AzAa - (A - l)c2) qA{AzX3 - (A - l)c2) ; ' 



det 17(t2, v) = A'^-^ (I) A{u, u), 



(53) 
(54) 



where 



A{u,v) = A 



AzAs - (A - l)c2 



zAa - c2 



AzA3ii2 Az2u2 
h2T + ^^(^i^3-y') 



(zAs - c2) 



(A - l)u 



r,2 



(a2c2 - a.c2), 



and 



^,(■5) = (AzAg - (A - l)c2) a, - (Azy - (A - l)a.c) Cj, 
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-j g= (zAi-a2)(zA3-c2)-(zy-a.c) 



(55) 
(56) 
(57) 



At this point we can calculate 



U=IJ, 



, which once completed will be substituted into equation (Eq). Firstly, 



we notice that the only variables which depend on v in the expressions for the I^^ are a and (3ki (defined in ([45D). 
Our first step therefore, is to simplify these quantities and then calculate their derivatives. 
Inserting equation (A42) into the expression for a and differentiating, we find 



dvidvj 



2q 



Z\d+1 92 



(zAa — c2) VA/ dvidv. 



(58) 



(From this point on we will use the notation g" to indicate that we are evaluating the differential at = /x.) 

To simplify (3-ki, we substitute for h and from equations ( |39| ) and ( A28| ) respectively; using equation (A40) we find 
that 
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2q 



ak 



zy — a.c 
zXn - c2 



We now note that, from the definition of b (equation ([A8|)) 



(zv - b.c) Q - (zAa - c^) 6z 



= 0, 



so the derivative of (5ki is given by 



dvidvj 



2qz 



(zAs — c^) dvidvj 



hkci. 



(59) 



(60) 



(61) 



the 



Differentiating equations (49)-(p2|) and using equations (|5J), ( p8| ) and (|6l|), we obtain the following expressions for 



^2 jii 



dvidvj 



dvidvj 



dvidvj 



dvidvj 



d7r(zA3-c2)[gA(l,0)]i/2 VA/ dv,d. 



d+2 ^2 



Z\rf+2 52^ 



^1 (1) 

tt{z\3 - c2) VA/ 
z'^SijhkCi V 



dvidvj 



rf«r!--\l,i)) 



A('i-2)/2[^A(l,5)]l/2' 



d7rA(zA3 — c2) dvidh 



-z^ hkCi 



7rA(zA3 — c2) dviduj 



1 J^^/(u^ 
^"[gA(ii,0)]i/2' 



(62) 
(63) 
(64) 



A('^-2)/2[gA(u,'C)]i/2 



(65) 



At this point we shall pause to give an overview of the simplification of equations (|62|)-(|65|). We first consider 
equations ( p^ ) and (^); on completing the contraction over the k and I indices, we obtain expressions in which the u 
integration may be completed exactly. We then complete the contractions over the i and j indices in the expression 



dvidvj 



E 



92 jy" 



,. duidvj 

fJ' n—l 



(66) 



and substitute this back into the expression for Ci{r, ti,t2) (equation |2^ ). 

Before we complete the contraction over the k and I indices in equations d64) and (pa), we will derive some results 
which will be used later. Using the definitions of h, A, k and q (equations (p9[), (p5[), (|56|) and ( p7| ) respectively), we 
find 



k({i).h= ^ A(u,^)-A 



AzAa - (A - l)c2 
zAa - c2 



(67) 



Using equation (|67| ) together with the definitions of Q ^{u,v), h and k (equations (|53|) , ( p9|) and ( ^6|) respectively) 
we also obtain: 



Xki{v) 



z{zX3 — c2)A(u, v) ' 

A 

z(AzA3-(A-l)c2) 



2^ k({i).CM2 

("^^"")^^'"A^A(^^^(") 



(68) 
(69) 
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Substituting equation (|68|) and ( p9[ ) into equations (64) and (|65|), we can therefore complete the contractions over the 
k and / indices to obtain: 



zk(0).c(5, 



dvidvj 



z V 



du 



A(^i, 0)3/2' 



(70) 



7rgi/2(zA3 — c^Y dvidvj 
Ak(^).c 



Az 



C^Cj 



dudv 



A('i-2)/2A(u,w)3/2 



3ki{v)kj{v) V? 



AzXs - (A - l)c2 gA(u, ^)(AzA3 - (A - l)c2) 



A(zA3 ~ c^)(Qfc.,-(i;) + c.,-fci(i;)) 
^(AzAa- (A - l)c2) 



(71) 



Equations (|7y) and (|7l| ) are now in a form in which we can complete the ii integration. We see that the only integrals 
required are A({t, C)^3/2 g^j^^j duu'^A{u,v)~^^'^; since A(u,{;) (which is defined by equation (|55|)) is of the 

form 7 + Su^, these integrals over u may easily be completed to give: 



du 



du 



(2A3-C2) 



[A(u,{i)]3/2 A(AzA3 - (A - l)c2)[A(l,u)]i/2' 



[A(u,{))]5/2 3A(AzA3 - (A - l)c2)[A(l,{))]3/2 



(72) 
(73) 



where 



gA(l, i)) = A2z2(AiA3 - t/2) - A(A - l)z(A3a^ - 2ya.c + Aic^) 
+ (A-l)2(a2c2-a.c2). 

We can therefore substitute equations (|7|) and ( [73| ) into the RHS of equations ( [70| ) and ( |7l| ) to obtain: 

k(0).c% a^?; 



dvidvj 



d7rA3(zA3 - c2)[gA(l,0)]i/2 dv.duj 
z d'^ V 



(74) 



(75) 



7r(zA3 — c2) dvidvj 

f k(^).c».-.i(l,f)) 
i (AzA3 - (A - l)c2) 



A''/2[gA(l,i))]l/2 

A(zA3 " c^)(fcj({i)cj- + Cjkjiv)) 
z(AzA3-(A-l)c2)2 



(76) 



We now substitute equations (|62|), (|63|), ( |75| ) and (|7^) into equation (|6^) and complete the contraction over the i 
and j indices. The details of this calculation are contained within appendix Following the contraction over the i 
and i indices, equation (pq) reduces to 



dvidvj 



A 

■KZ 



poo 6 



(77) 



where the terms T„ are defined by the equations ( |C28| )-( C33 ). Substituting equations ( |20| ) and (|77| ) into equation 
(p8|), where without loss of generality we can set x = 0, we obtain the following expression for Ci{r, <i, t2): 



dr 



nz{4:TTD{t2 - t))'*/2 



AD{t2 



(78) 



The spatial integral over fj, may be completed by transforming to spherical polars and choosing a change of variables 
which allows the integral to be completed by steepest descents. In spherical polars the spatial integral may be written 
as 
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/ 



2(n) 



r(^) 



dip (sin V') 



rf-2 



(79) 



We notice that each expression for T„ (equations ( C28 )-( p33|) ) contains the factor A"''/^ = (1 + 2\v j z)~'^l'^\ we 
therefore make the substitution w = vXd/z, so that in the large-d hmit, this factor reduces to an exponential in 
w. We also make the additional substitution r = t2{l — x/d), and rescale ii to fi — {xt2)~^^^ fJ.- The form of the 
r substitution is chosen in the expectation that the contribution to the integral from later times will dominate the 
result. Equation (78) then reduces to 



Ci{r, ti,t2) 



4*2 



7r3/2d2r(£tJ.) \ADJ 

I 

dx 



, \ d/2 6 

—) y 

AD / ^ 



dwi f 



dip 



A Jo (sintpy 



■exp[-dg{jj,ip)] T„, 



(80) 



where 



ln(/i sinip). 



(81) 



Once the integral is in this form, the fi, and ip integrals can be completed in the large-d limit by steepest descents. 
Applying this limit does not represent a serious limitation on the calculation, since this is the regime of most interest, 
as explained in the introduction. In the limit d oo, provided the large-d behavior of the terms Ti is controlled (as 
demonstrated in the first section of appendix Q), the value of the integral is dominated by the contribution from the 
neighbourhood around the minima of the function g{p,ip). 

Within the range of integration, the function g(fi, ip) has a global minimum at /i = \/'2D, ip — 7r/2. We can expand 
the integrand to first order about this minimum, reducing the jl and ip integrations in equation (pQ) into Gaussians, 
and giving 



Ci(r,ii,t2) = 



4*2 exp(-d/2) 
(2DY/^n^l^d^T{^) 



d\ 
2 



)d/2 6 noo 
y dx dw T„| y2jj ,^^^/2 

d(V-7r/2)2^ 



dpL exp 




dip exp 



(82) 



Since we are only interested in the first order correction to the correlation function, we use Stirling's formula to 
evaluate the leading order contribution from the Gamma function, which is given by: 



1 



2(2^)1/2 fd 



d/2 



exp(-d/2). 



On substituting equation (R3|) back into equation (p2[) and completing the Gaussian integrals, we find 



Ci(r,ii,<2) = 



2,t ^ 



(83) 



(84) 



Now all that remains is to calculate the leading order contribution from the integration over x and w of the terms 
T„, evaluated at /i = ^/2D and ijj = 7^/2. Each of the terms may be expanded as a power-series in l/d, and in section 
2 of appendix ^ each term is evaluated to 0(1), giving: 



X 1 



7(*2 — ti) exp(— a;/2) w exp(— tD) 

%D^h+t2){l^^lW^ 

7(*2 — ^1)2; exp(— a;/2) w exp(— w) 
7(^2 — ti)x(x — 2) exp(— a;/2) exp(— ?«) 

%Dt2{u + t2){i-i^y/^ - 



2*27 r 



2„2 



4*1*27^ 
(*l+t2)V ' ^(^1+^2)2 



(85) 
(86) 
(87) 
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-'yr'^x exp{—x/2) exp{—w) 
-7 exp(~a;/2) exp(— it) 



"t2 - h t2r^ 



X 1 



2*27 I-' 



2„2 



4*1*27 

(*l+t2)2; * i?(*l+t2)^ 

4*2 



*1 +*2 I?!*! +*2)2 2 V (*1 +*2) 



We now complete the x and t« integrations of the expressions given in equations (|85|)-(90), finding: 



oc roo 



'Q JO 



dx dw Ti 
dxdw T2 



7(*2-*i) 



Jo 



00 /"OO 

dx dw T-i 



dx dw T4 



2L'*2(*l +*2)(1 - 


-y2)3/2 


7(^2 -*i) 




£'*2(*1 +*2)(1 - 


^2)1/2' 


7(^2 -<i) 




£'*2(*1 +*2)(1 - 


y)l/2' 


9 

— 7r 





4*1*27^ 



2„2 



*27 r 



2i?2(*i +*2)2(1~ 72)1/2' 



dxdw Tk = 



-7 



2i:'*2(i - 72)3/2 

<2-<l 



1 - 



4*2 



2- 



8*1*27^ 



2„2 



*27 



*2r2 



(^1+^2)2; V (^1+^2)2 i?(*l+<2)2 

4*1*272 ^ *272r2 



*l+*2 £>(*!+ *2)2 



(*l+*2)2 i?(*l+*2)^ 



dx dw Tr — 0. 



On substituting equations (91 )-(|96[) into equation (B4), we obtain: 



C'i(r,*i,*2) 



27 



£l7rd2(l -^2)3/2 



4*2 



4*1*27^ 

(*i + *2)2; V(ii + ^2) 



2(*2-*i)(l-7') , 2*2(*2 -*i)72r2 



(*1 + *2) 



i?(*l+*2)4 2D-\ti+t2Y 



(88) 



(89) 
(90) 

(91) 
(92) 
(93) 
(94) 



(95) 
(96) 



(97) 



Finally, on substituting equation ( |97D back into equation (|25D, we obtain the complete expression for the first order 
correction to the correlation function: 



Ci(r,*i,*2) = 

' *! -*2 
*1 +*2 



47 



1377^2(1 -^2)3/2 
2 / , 2 



1 



4*1*27^ 



9 9 



(*l+*2)2 D{ti+h)J 4D^h+t2) 



(*2+<2)72r4 



(98) 



We now express our final result for the first order correction to the correlation function in the same form as the 
exact OJK result (equation (jl^)). In this form the correlation function is given by 



C(r,*i,*2) = -sin-i[7(l + AF(r,*i,<2))], 

IT 



(99) 



where 
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{tl+t2f 



2 ? 



(100) 



We can clearly see from this expression that the first order correction to the correlation function is 0(1 /d?), which 
lends weight to the assertion that the OJK result becomes exact in an infinite-dimensional system. 



A. Two special cases 

We now evaluate this result in two special cases: at equal times, and when the times are widely separated, ti 3> ^2- 
At equal times, the first order correction term is given by 



where 



Ci(r,i) 



1 / r 
t9 



exp(-3a;V8) 



(101) 



(l-exp(-a;2/4))3/2' 

This correction term clearly exhibits the expected scaling, L ^ i^/^; the scaling function g{x) is shown in figure |^. 



(102) 




FIG. 1. Scaling function for the first order correction to the correlation function. 

Figure || shows a comparison of the zero order OJK result (equation (|l5|)) with the perturbed result for d ^ 2,3 and 
4 at equal times; the functions have all been scaled so that they have the same gradient at the origin. Although the 
result is only valid for large d, the figure clearly demonstrates that the perturbation will have the effect of lowering 
the exact OJK result. This is discussed further in section [V, 
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X! 





FIG. 2. Comparison of the OJK and perturbed correlation functions, shown for d = 2, 3, 4. The solid line is the OJK result, 
which is independent of dimension. As the dimension increases, the size of the perturbation decreases. 



If we now expand equation (101) about r = 0, we see that for small r, 

^■C' " - 5^ (ZJ^VI (' - iSl + SlS^ + (103) 

and hence the result obeys Porod's Law and the Tomita sum rule ||l^,|l^ ■ 

We now consider the case where the times are widely separated. If ti 3> ^2, then equation (|98| ) reduces to 



Comparing this to the conventional scaling form, 

Ci(r,ii,f2) - (^^^ ''(zt) *i»*2, (105) 

where Li and L2 are the characteristic lengths for the system at times ti and ^2 respectively, we find that A = (i/2, 
as expected. We also note that the leading order correction to the correlation function in this limit has exactly the 
same form as the zero order correlation function; from equation ( p^ we have 

2/4^2^''/' 



Co(r,tr,t2) = -(-l expl--^). (106) 



IV. DISCUSSION 



In this paper I have described a perturbation theory applicable to OJK theory in a system with an infinite number of 
dimensions and no noise. The leading order correction term to the two-time correlation function has been calculated 
by treating the non-linear terms in the OJK auxiliary-field evolution equation as a perturbation to the linearized 
equation. There are two separate motivations for this study. Firstly, since in conventional OJK theory, the evolution 
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equation is linear and Gaussian initial conditions are imposed, it follows that the distribution of the auxiliary field must 
be Gaussian at all times. This assumption, present in many of the approximate theories, has been critically assessed 
by Yeung et. al. they show that the auxiliary field distribution measured directly from numerical simulations is 
not exactly Gaussian. Hence, the advantage of the approach outlined above is that the auxiliary field distribution 
may take any form. 

Secondly, it has been proposed by several authors that OJK theory is exact in an infinite number of 

dimensions. Determining the dimensional dependence of the correction term allows this hypothesis to be examined. 

The main result of this paper is that the first order correction term to the correlation function is 0{1/(P), which 
gives further confidence to the assertion that OJK theory becomes exact in the large-d limit. This may be compared 
to the calculation of Liu and Mazenko |]l3| in which they made a perturbative expansion near d = 1 and d = oo 
about the approximate theory developed by Mazenko [p^-pO|. At large d, the first order correction term in their 
expansion was found to be 0(l/d). The relative sizes of these corrections are consistent with numerical simulations 
1^ which show that the OJK result provides a marginally more accurate prediction of simulation data for the pair 
correlation function than the more sophisticated approach of Mazenko. However both results do suggest that the 
results of OJK are asymptotically exact in infinite-dimensional systems. We also note that in a recent work, Mazenko 
pO| has developed a perturbation expansion in the cumulants, which at zeroth order recovers the OJK result. At 
second order, however, significant deviations from the OJK theory are obtained at large d, in contrast to the 0{l/d^) 
correction obtained here. 

The two-time correction term to the correlation function is given by equation (p8|). In the limit ti 3> t2 (on 
comparing equations (104) and ( |106| )), we find that the correction term has exactly the same form as the zero order 
result, and the expected OJK result, A = d/2, is recovered. 

If we examine the leading order correction at equal times, given by equation (101), the first observation is that this 



function is odd in r, as is the zero order term. Porod's law [ [14[ is therefore obeyed; the 0{r) term in the expansion 
(equation ( |103[ )) modifies the amplitude of the large-fc tail, which is proportional to the density of defects ||l|,| 
absence of an term ensures that the correction term satisfies the Tomita sum rule Hjie), dk{k'^ 
where A is the amplitude of the Porod tail. 

We also note that the perturbation has a negligible effect on the large-distance behavior of the correlation function, 
see figurejl|. 

Figure ||demonstrates how the perturbation modifies the OJK result, plotting Co -t- Ci for d = 2, 3, 4. Although this 
calculation is only valid at large d, this graph demonstrates that the correction term will have the effect of lowering 
the OJK result. This is exactly as expected since Humayun and Bray |^ showed by comparing simulation results 
with OJK theory, that while the OJK result provides an accurate prediction for initial conditions with short-range 
correlations, the theoretical result is slightly higher than the simulation data. However, it has also been demonstrated 
that when long-range correlation are present in the initial conditions, the OJK results are no longer satisfactory 
This suggests that a possible extension to this work could be to consider the effects of long-range initial conditions 
on the calculation. 

The main limitation of this calculation is that to retrieve the full evolution equation for the auxiliary field (equation 
(^), we need to set the perturbation parameter, A to 1. This means that the calculation of the entire correction term 
requires a sum over all orders in A. However, due to the complexity of the present calculation, I have been unable 
to evaluate the correction terms at higher orders in A. In principle the sum of the higher order terms could alter the 
d-dependence of the correction term, but the result remains a strong indication that OJK theory becomes exact in 
an infinite-dimensional system. 

Finally, since we noted that the distribution of the auxiliary field is non-Gaussian, it is of interest to consider the 
exact form of this distribution. This can be calculated as follows; the distribution can be written in the following form, 
P{x) = {S{x — m{r,t))), this may be expanded in A using m(r, t) = mo(r, t) + Xmi{r, t) + 0{X'^), to give P{x) = {5{x — 
mo(r,t))) — X{mi{r,t)-j^5{x — mo(r, t))). The first term on the RHS reduces to the expected Gaussian distribution; 
the first order term in A can be evaluated using a similar method to the correlation function calculation, i.e. by 
inserting the expression for mi(r, t), multiplying by the relevant probability distribution function, and completing the 
integrals. However we will leave this question to future work. 
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APPENDIX A: EVALUATION OF THE JOINT PROBABILITY DISTRIBUTION 



In this appendix the joint probability distribution is explicitly calculated. It is defined by 



P(m) = 



1 



(27r)^(detA-i)^ 



exp 



" 2 



(Al) 



where A-^^ — {rhifhj) and the vector m is defined in equation (^2|). The first step is to calculate all the correlators 
which define the elements of the matrix A^^; we can then find the inverse of this matrix {Aij) and the determinant, 
det^-i. 



1. Calculation of the Correlators 



In this section we evaluate all the elements of the matrix A^^^ = {rriirhj): for this we require the following correlators: 
(TOo(x,i)mo(x',i')): {mo{x,t) ^^mo{x.' ,t')) , and (g|-mo(x,t) -£^mo{x,t)). 

Substituting equation {Uw into equation (18) gives an expression for the auxiliary field too(x, t): 



mo(x,t) 



1 



d^/i mo(/i, 0) exp 



(x-m) 



2 1 



and substituting this expression into the correlator (too(x, t) mo(x', t')) gives 

1 



(A2) 



(mo(x, t) mo(x ,i )) 



exp 



(x - (x' - rif 



(A3) 
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where we have aheady appHed the conventional Gaussian initial conditions, 

(too(x + r, 0) mo(x, 0)) = A<5(r). 
This integral may be evaluated by completing the square in the exponent and making a change of variables /x' 
— {t'x + <x') /{t + t'), leaving a simple Gaussian integral which, once completed, gives 



(mo(x, i)mo(x',t')) = 



A 



exp 



(x - xQ^ 



The remaining correlators are easily obtained by differentiating equation (A4); we obtain: 

(TOo(x,i) ^mo(x',i')> - 2D{t + i) ("^o(x,t)TOo(x',t')), 



Substituting equations ( |A4| ), 
(Ar.^ = {rhirhj)), we find 



and (A6) into the definition oi A ^ 



(A4) 

(A5) 
(A6) 



= (-) 

z 



/ 


1 











: \ 










a 


b 


c 









1 
















zXi 


zw 


zy 










zw 


z\2 


zv 


V 








zy 


zv 


z\J 



(A7) 



where: 



A 



z = ADt, 
A = ■ 



37r_Dr)<i/2 ' 



A2 - 1, 
h = {v — /x) exp 



SDt 



w={^) exp 

« = (ti^) ' exp 
a=(x + r-/x) (^) 



4D(ti+r) 



4n(«i+t2) 



4D{t2+T) 

{d + 2) 



exp 



(x+r-/X)^ 

4D(ti+r) 



C = 



exp 



4-D(t2+r) 



(A8) 



2. Calculation of A 



The elements of the matrix A are calculated by constructing the adjoint and determinant of the inverse, since 
A = Adj(A~^)/det(A~^). Let the elements of the adjoint be defined by: 



Adj(A-i) = 



V. 



Oij 




V 


C 


e ■ 


■ p 


s 


t 




s 


q 


u 


e ■ 


■ t 


u 


r) 



(A9) 



The elements p, g, r, ■ 
matrix of the form 



■ , u can then be calculated directly using a formula for the determinant of a (d + 2) x (d + 2) 
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1 1 



B 



a' b' 



1 : : 

••• c'^ ••• p' q' 



(AlO) 



r' s' I 
which is given by 

detB = (p' - a'.c')(s' - b'.d') - ((?' - b'.c')(r' - a'.d'). 
We derive this result by contracting the free indices in the equation 

det-B = ^ eij,,,xyzBiiB2j ■ ■ ■ Bd+2z, 



(All) 



(A12) 



i]...xyz 



where €ij...xyz is a {d+2) anisotropic tensor, which takes the vahie 1 if (ij . . . xyz) is an even permutation of (1 2 3 . . . {d-\- 
2)), —1 for an odd permutation, and zero otherwise. 

First consider the sum over the index i : Bi j is only non-zero ifi = l,rf+lor(i + 2, [Bi i = 1, Bid+i = a'l and 
Bi d+2 = b'l), which implies that equation (|A12|) may be written as 

(A13) 



deti? — ^ ilj...xyzB2j ■ ■ ■ Bd+2z + Di, 
jk...xyz 



where 



k...xyz^\B2 j ■ • ■ B(iJ^2 z ^d^2jk...xyz^\^2 j • • • ; 

jk...xyz 



(A14) 



We now complete the sum over j in equation (A14) by noting that the first term on the RHS of this expression will 
only be non-zero if j = 2 or j = d -|- 2, whereas the second term will only make a non-zero contribution if j = 2 or 
j = rf + 1, giving 



^d-\-ld-^2k...xyz(^l^2^^ k • ■ • ^d-\-2 z ~l~ ^rf+2d+lA;...a:y£^l'^2-^3 k • ■ • ; 

k...xyz 



-\-^d'^12k...xyzO'iB3k ■ • ■ Bd+2z + ^d^22k...xyzb'iBzk • • • Bd 



+2 z 



(A15) 



We can now complete the sum over all the remaining indices for the first two terms on the RHS of equation ( [A15| ), 
since for a non-zero contribution we must have fc = 3, / = 4, . . . , a; = d and eit her y — 1 and z = 2, ov y — 2 and 
z = 1. After completing the sum over k on the remaining RHS terms, equation ( A15 ) gives 

+ ^ ( f-d+12d+2l...xyza''lb'j,Bii . . . Bd+2z + £d+23d+ll...xyzb'ia'^B4i . . . Bd+2z 
I. ..xyz 



+ ^d+123l...xyza.'iB4i . . . Bd+2z + ^d+223l...xyzb'iBii . . . B, 



d+2 z 



(A16) 



We can use this method repeatedly to sum over all the free indices in equation ( A16 ), eventually giving 



—sa'^c'i + qa[d'i ~ pb[d'i + rb'ic'i. 



(A17) 



Substituting equation (A17) back into equation (A13), we get: 
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detB = ^ eij...xy2-B2 j ■■■ Bd+2 z - sa'ic'i + qo!^(i^ ~ pb[d[ + rb[c[ 

jk...xyz 
d 

We can now repeat this entire procedure to complete the sums over the indices j . . .x, obtaining 
detS = ei2. ..dyzBd+iyBd+2z - sa'.c' + ga'.d' -pb'.d' + rb'.c' 



i—1 j>i 



Two of these terms simphfy further; we find: 

y^^ei2. ..dyzBd+iyBd+2z = ps - qr, 



(A18) 



(A19) 



^ J2 {o^Ab'A + Kd^'^A ~ <dfJ^A " K^o/^d'^) = (a'.c')(b'.d') - (a'.d')(b'.c'), 

i=l j>i 

and therefore equation ( A19| ) finally reduces to 

detS = {p' - a'.c')(s' - b'.d') - {q' - b'.c')(/ - a'.d'). 



(A20) 



Having calculated the determinant of B, we can use this result to evaluate the elements p, g, r, • • • , m of the adjoint 
matrix (equation (|A9|) ) . We obtain 



X 

z 

A 

z 

A 

z 

A 

z 

X\ 

z ^ 

A 

z 



d+2 



d+2 



d+2 



d+2 



d+2 



d+2 



[(2A2 






— {zv — 


h.cf] , 


[(zAi 


-Si^KzXs 




- {zy - 


a.c)^] , 


[(zAi 


-a^)izX2 




— (zw - 


- a.b)^] , 


[{zv - 


'b.c)(zy- 


a.c) - 


- [zw — 


a.b)(zA3-c2)] 


[{zw - 


- aL.h){zv - 


b.c) 


- {ZX2 ^ 


-b2)(zy-a.c)] 


[{zw 


- a-h){zy - 


a.c) 




- s?){zv - b.c) 



(A21) 
(A22) 
(A23) 
(A24) 
(A25) 
(A26) 



The remaining elements of the adjoint are evaluated directly by a component-wise expansion of the equation 
Adj(A^^) = dci{A^^) I, giving the following expressions: 



I = -(pa + sb + te), 
T] = — (sa + qh + uc), 
C = — (ta + uh + re), 



%j = -(detA ^)6tj - flj^j - b,T]j - CiQ. 



(A27) 
(A28) 
(A29) 

(A30) 



Now that we have evaluated all the components of the adjoint, we can evalu ate the expression fo r th e joint probability 
distribution, -P(m), in terms of these variables. Substituting equation ( [A9| ) into equation (Al) and using A = 
Adj(yl-i)/det(A-i); we find 
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P(m) 



i^(m) 



2detA- 



where 



F(m) = ( m'kOkim'i + 2|.m'm(l) + 2T7.m'm(2) + 2C.m'm(3) + pm(l)^ + qm{2f 
+rm(3)^ + 2sr7i(l)m(2) + 2tm(l)rn(3) + 2MTO(2)m(3) 



(A31) 



(A32) 



3. Identities Relating the Components of the Adjoint oi A ^ to the Determinant of ^ ^ 

In the calculation described in this paper, there are several expressions which are frequently used to simplify the 
algebra. These are all derived directly from the component-wise expansion of the equation A^^ Adj(A~^) = detj4~^ /, 
and are listed here for ease of reference. 



p{z\i — s?) + s(zw — a.b) + t{zy — a.c) = — detA 

A 

s(zAi — a^) + q{z'uj — a.b) + u{zy — a.c) = 0, 
t{z\i — s?) + u{zw — a.b) + r{zy — a.c) = 0, 
p{zw — a.b) + s{z\2 — b^) + t{zv — b.c) = 0, 

s{zw — a.b) + q{z\2 — h^) + u{zv — b.c) = -^dctA""^, 

A 

t{zw — a.b) + u{z\2 — b^) + r{zv — b.c) = 0, 
p{zy — a.c) + s{zv — b.c) + t{z\^ — c^) = 0, 
s{zy — a.c) + q{zv — b.c) + u{z\i, — c^) = 0, 

t{zy — a.c) + u{zv — b.c) + r-(zA3 — c^) — ^ActA^^. 

A 



(A33) 

(A34) 
(A35) 
(A36) 

(A37) 

(A38) 
(A39) 
(A40) 

(A41) 



I n ad dition, by using the definitions of p, g, • • • , u (from equations ( A21 )-( A26 )) and equations (A32), (A34), ( [A37| ), 
(|A39|) and (|A40|), we can see that 



{pq-s^)^ 0) ^ detA-i(zA3-c2), 



{qr - u^) 
(us — qt) - 



d+l 



d+l 



I detA ^zy-a.c). 



(A42) 
(A43) 
(A44) 



APPENDIX B: CALCULATION OF THE INVERSE AND DETERMINANT OF Q. 



In this appendix, the expressions for the determinant and inverse of the matrix fi, which are required for the 
completion of the integrals in equations (E9|)-(p3), are calculated. The matrix Vl{u,v) is defined by 



1 



AetA 



-1 \ 



^ + '^K ^hAu^ - 1) ) + 2^5,, 

q q I - 



(Bl) 



where h, r] and 9 are defined by e quati on s (|39| ), ( A27) a nd ( A30| ) respectively. Substituting for h, 77 and 6 into 
equation (Bl), and using equations ( A22 ), ( A29| ) and ( A42 )-( A44 ), we find the expression for Q{u,v) reduces to 



A Sij + fi aiUj + ^ [ttiCj + CiGj) + i> CiCj 



(B2) 
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where: 

A = l + 2Au/z, (B3) 

/i=^(zA3-c2), (B4) 

e = -^(z2/-a.c), (B5) 



(i+2 



q = (zAi - a^)(zA3 - c^) - (zy - a.c)^ (B7) 



1. Determinant of Q.{u,v) 

The determinant of the matrix f2(M,w) is calculated by evaluating the product of its d eigenvalues. Any vector 
orthogonal to both a and c will have eigenvalue zK/X^ so we need only calculate the two remaining eigenvalues, which 
are associated with the eigenvectors which lie in the plane spanned by a and c. The eigenvalue equation for these two 
may be written as 



(B8) 



We can substitute for ^(u, v) from equation (Bl) and complete the sum over the index j; by equating the coefficients 
of ai and Ci we then obtain two simultaneous equations for [3 and 7, given by 

A + /ia^ + ^a.c + 7(/ia.c + ^c^) 



Pi 



^a + i>a.c + 7(A + ^a.c + i^c 



(B9) 
(BIO) 



On eliminating 7 between these two equations, we have a quadratic equation in from which we extract the product 
of the two roots of the quadratic, P±. This product is given by 



where 



/3+/3_ = (^-j A{u,v), 



A{u,v) = (A + ^a.c + i>c^)(A + /ia^ + ^a.c) - (^a^ + i>a.c)(/ia.c + ^c^ 



(Bll) 



(B12) 



To simplify this expression we expand the brackets on the RHS and use equations (B4)-(B7); after some algebra, 
equation ( |B12 ) reduces to 



A{u,v) = A 



AzAg - (A - ly 
zh - c2 

(A-l)fi2 2 



+ ^^(AiA3-y") 



(zAs - c2) 



q 



(a^c^ - a.c^). 



(B13) 



We recall that the remaining eigenvectors all have the eigenvalue zA/A, and therefore the product of all the eigenvalues, 
which is equal to the determinant, is given by 



detn = A 



d-2 



A{u,v). 



(B14) 



20 



2. Inverse of fl{u, v) 



We begin to find the inverse of n{u, v) by defining the variables A, B, C and V by the equation 

fl^-^ — —{ASij + BaiGj + C {ttiCj + Ciflj) + V CiCj). 



(B15) 



This must satisfy the identity flij fljf^ — Sik- Hence, by expanding this equation and equating coefficients, we obtain 
the set of simultaneous equations given below: 



0(1) AA = 1, 

0{a,ak) Afl + B{A + jla? + ^a.c) + C(/ia.c + ^c^) = 0, 

0(a,cfe) AS, + C(A + fia? + Ca.c) + P(^a.c + ^c^) = 0, 

0(c,afc) AS, + B{ia? + v&.c) + C(A + |a.c + vc^) = 0, 

0{c^Ck) Av + C{is? + vsi.c) + V{A + |a.c + i>c^) = 0. 



(B16) 
(B17) 
(B18) 
(B19) 
(B20) 



Since we have one more equation than we require to determine t he so lu tions for A, B, C and V, we discard one 
equation and check f or co nsistency later. Solving equations (B16), ( B17 ), ( B18 ) and ( B20 ) simultaneously, and using 
equations (B7) and (B13) to simplify the results, we find, after some algebra, that: 



AqA 



A = 
B 

C 

V 



AqA 
AqA 



(AzA3-(A-l)c2), 



(Azy- (A - l)a.c). 



(AzAi - (A - l)a2) + 



ju' 1) 
A(zA3-c2)- 



(B21) 
(B22) 

(B23) 

(B24) 



Substituting these results back into equation ( ]B15| ), we see that the inverse of Q{u,v) is given by 



A{U^ - 1) /A 1N N 

+ -rT;{Azy — (A — l)a.cj(ajCj + Ciaj) 



qA 



A(zA3 - c2) ' qA 
((AzAs - (A - l)c^)aiaj + (AzAi - (A - l)a^)ciC-,) 



(B25) 



This expression may be simplified considerably by the introduction of a new variable. If we define k by the equation 

h = {AzXs - (A - l)c2) a, - {Azy - (A - l)a.c) c^, (B26) 
then equation ( |B25| ) may be rewritten as 



(AzAa - (A - l)c2) ^A(AzA3 - (A - l)c2) 



(B27) 



APPENDIX C: CONTRACTION OVER THE / AND J INDICES 



In this appendix we calculate a non-indexed expression for 



dvidvj 



M ^1 ^^^^^^ 



(CI) 



where: 



21 



dvidvj 



dvidvj 



52 



dvidvj 



dvidvj 



H JO 
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7r(2;A3 — c2) dvidvj 

k(0).c% 

d^A3(zA3 - c2)[qA(l,0)]i/2 dv.dv, 

dv 



A(''-2)/2[gA(l,5)]l/2' 



(C2) 
(C3) 
(C4) 



(0A3-c2)yo A'i/2[^A(i,w)]i/2 ajy-ai.^. 



k{f}).cr!,/(l,{)) A(zA3-c2)(fc,({i)cj+c,%(i)))' 



(AzAg - (A - l)c2) 



z(AzA3-(A-l)c2)2 



(C5) 



and s = (z/A)''+2s^ g being defined by equation ( A24| ). 

To obtain such an expression, we first evaluate the sum of the four expressions given by equations ( |C2|) -(C5), and 
then contract this expression over the free indices i and j. 



To simpfify these expressions we need to calculate the second derivative of s. Using equations (A24) and (|60|), we 
find that 



925 



dvidvj 



= z (2?/- a.c) 



dvidvj 



{Z\z - c2) 



dw 



dvidvj 



(C6) 



(C2)-(|C5|), we may combine the resulting expressions to obtain: 



After substituting for k and the second derivative of s (from equations ( pq ) and (C6) respectively), into equations 

I , (C7) 



dvidvj 



92 {i^'+ri) 



7rdA3(gA(l,0))i/2 I ^ diy.diy. 



dvidvj 



A3 



92 ui 



dvidvj 



fJ'/ 



dv 



7rA'*/2(gA(l,i}))i/2 



A {k{v)iCj + Cik{v)j) 92 V 



(A2A3 - (A - l)c2)2 dv^dvj 



where 



R 



A3 / Kzy - (A - l)a.c 



y \Az\3 - (A - l)c2 



( yR 






Jl 


1^A3 


dvidvj 






= 1 + 


(A-l)(yc2-A3a.c) 




y(A2A3 


- (A - l)c2) 





(C8) 



(C9) 



Before attempting to calculate the sum in equation (p4|), we manipulate equation (C8) so that part of the v integral 
can be completed exactly. This will simplify the algebra greatly, since the exactly integrable term in equation (CS) 
will cancel the contribution to the sum from equation (C7). 

We define a new variable "^ij , by the equation 



A 



Az 



(CIO) 



on comparing this with equation ([B25| ), we find: 

'I'jj — {AzXz — (A — l)c2) QiGj — [Azy — (A — l)a.c) {aiCj + ajCi) 
+ (AzAi - (A - l)a2) c,Cj. 



(Cll) 



Substituting for R and O- ^ from equations (|C9|) and (CIO), equation (|C^) can be written as: 
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= s-- 



dv 



TT Jo A'i/2(gA(l,{)))l/2 



{k{v)iCj + Cik{v)j) d^v 



(AzAa - (A - l)c2)2 dv.dvj 



(A-l)(yc2-A3a.c),5,,- d^v 
^ A3(AzA3-(A-l)c2) du.dy, 







92 u; 


.)1 


dvidvj 




dvidvj 





where 



5 



A(5y / y 52^ 



A'i/2(gA(l,w))i/2 



TT \ A3 dVidVj 

We now define a second new variable, Q, by the equation 

A(l,i}) = A2A(l,0)(l-Q(i))), 

and, using equation (Q), we find that 

A(A - l)z{Xj,a? - 2ya.c + Aic2) - (A - I)2(a2c2 - a.c2) 



A2z2(AiA3-y2) 



(C12) 



(C13) 



(C14) 



(C15) 



On substituting for A(l,-D) from equation ( pl4 ) into equation (C13), and expanding the expression (1 — Q) ^/2 
using the binomial theorem, we find that 



y d^v 



7r(gA(l,0))i/2 I A3 dv.dvj 



^ dvidvj 



My 



dv 



\{d+2)/2 



(CI6) 



Hence, using A ('^+2)/2 — (Ad), and inserting equation (C7) in the expression for 5*, we find that equation (|C16| ) 
reduces to: 



S = -- 



XSij 


y d^v 




32 w 




dv 


7r(^A(l,0))i/2 


A3 dvidvj 






.) 





E 



dvidvj 



(C17) 



Next, substituting equation (C17) into equation (|C12| ) and using equation (CI), we find that 



dvidvj 



dv 



7r(9A(l,0))i/2 7„ A('^+2)/2 Vm^ 



A 



dv 



7r(gA(l,0))i/2 7„ A(''+2)/2(l-Q)i/2 



(A-l)(yc2-A3a.c)% d^v 



A3(AzA3-(A-l)c2) dvidi^j 



(fc(u)jCj- + Cifc(w)j) d'^v 



(AzAs- (A - l)c2)2 ^l^^^L'J 



A2gA(l,0)(l-Q) I A3 diy^di 











dvidvj 




dvidvj 





(C18) 



Before we start the contraction of equation (C18) over the indices i and j it is convenient to derive some useful 
results. 

First, we calculate the derivatives of v and w (from equations ( [A8| ) ) , these are given by: 

(C19) 
(C20) 



92 tl 




CiCj 


vStj 




dvidvj 




z^v 


2D{t2+T)_ 








aiQj 






duidvj 






2D{ti+T) 
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(From this point on, we will use the notation; and t o rep resent v and w evaluated f = fi). 

Using these equations together with equations ( ^6|) and ( Cll ), we evaluate the following expressions: 



9^ V 




if 


dvidvj 


(J- 




d"^ w 




if 


dvidvj 


(J- 





5i. 

{k{v)iCj + Cik{v)j) 



Vf,d 



z^v^ 2D{t2 +t)J' 



dvidvj 



2zA(A3a.c - yc^ 



Z-^V^ 2D{t2+T)J' 



5ij = Az(A3a2 - 2ya.c + AiC-") - 2(A - l)(a^c^ - a.c^), 



C V \ 

(AzA3-(A-l)c2) 



(a^c^-a.c^), 



2D{ti + r) 



'A" 



(AzAi - (A - l)a2) 



z'^w^ 



(a^c^ - a.c^ 



Final ly, we can compl ete t he contraction over i and j by substituting equations ( p21 )-( C2^ ), together with 
( IcTsI) ) into equation ( |C18| ) to give: 



(C21) 

(C22) 

(C23) 
(C24) 

(C25) 

(C26) 
^) and 



where: 



duidvj 



= / dvS^Tn 



Ti 



^^-^ [I) 

(A1A3 - j/2)l/2A(<i+2)/2 
,2 



Wf^d 



A3 V^'^M 2D{t2 + T)) \z^w^ 2D{ti+T) 
(A-l)(;;c2-A3a.c) 

A(rf+2)/2^3(A^A3 - y2)l/2(l - 0)l/2(AzA3 - (A - l)c2) 

Vnd 



Z^V^ 2D{t2+T)J' 



2z(A3a.c - yc^) 



Arf/2(AiA3 - 2/2)1/2(1 _ Q)l/2(AzA3 - (A - l)c2)2 
,2 



Z^Vf, 2D{t2+T)J' 

z^^{a?c^ — a.c2) 

A(rf+6)/2(A^A3 - 7/2)3/2(1 „ Q)3/2 

/ (AzT/- (A - l)a.c) (AzAi - (A - l)a2) 
-(A3a2 - 2ya.c + Aic2) 

A(rf+4)/2z(AiA3 - y2)3/2(i _ g)3/2 

(y]^ 

I A3 \Z^V^ 2D{t2+T] 



z^Wf, 2D{ti+T] 



(C27) 



(C28) 



(C29) 



(C30) 



(C31) 



(C32) 
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2(A-l)(aV-a.c^) 

6 A(^+6)/2z2(AiA3 -?;2)3/2(l _Q)3/2 

? \ \ 

(C33) 




z'^Vf, 2D{t2+T)J \z^Wf, 2D{ti+T 



APPENDIX D: 



1. The Large-d Behavior of the r„ 



In t his se ction we will demonstrate that the large-d behavior of the terms, Ti, which are defined by equations 
( |C28D -( C35 ), is controlled, and hence that we are justified in applying the method of steepest descents to the jj, and 
■0 integrals in equation (pO|). 

To simplify the algebra within this section we define several new variables: 



d/i 



/2 = (1 - x/2d), 



E = f2 ^ exp 



exp 



(r - {xt2)^lJ-f 



AD{ti+t2) ^D{ti+t2)h J ■ 
Therefore, the set of variables given by equations (|A8|), with which the terms Ti are defined, may be rewritten as: 



= 7/3/1"" (If) 



z^ADt2h, 

d/2 



d/4 



'^3 — 7i , 
V^^f2ft^'E, 



2t27 j(d+2)/2 /t2\ 
ti+t2-ll \ti) 



d/i 



(Dl) 



r - [xt2)^ii 



c^-{xt2)ifr'^^'f.E. 



To investigate the large-d behavior of the r„, we substitute these variables back into equations ( p28 )-( C3^ ) to 
obtain: 



27:»A(rf+2)/2(l_ ^2)1/2 



E / xfi^ 



Es Mr - ixt2)^fi) 



ir.^2 



2t2 \ 4Df2 



h+t2 \ 2Dh{ti+t2) 



(D2) 



To 



-7(A-l)/('^+2)/^£;2 



xn 



T3 



ADt2{tl + t2)A(<^+2)/2(l - 72)1/2(1 _ Q)l/2 yADf2 

{x(i^[2t2Es - iti+t2)E] - 2{xt2)^ fi.rEs) 
ADA - (A - l)//'^^^^^xfi^E-^ 

Ad/H2{h + t2){l - 72)1/2(1 - Q)i/2 yADf2 J 
{xfi^[2t2E, - {ti+t2)E] - 2{xt2)^ (i.vE,) 



(D3) 



(4DA-(A-1)/} 



(<i+2)/2 



xfi^E'^f 



(D4) 
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where: 



(d+2)/2 



EEs (A^r2 - fi.v 



16i:>4A('^+6)/2(ti + t2Y{l - 72)3/2(1 _ g)3/2 
'l^E: 



12 



DK - ^ J/.^ . . ^(a;t2A - {xt2)^ (i.r) 



jE 



DA 



2{h + t2) 



.^(l_^2)^3/2^(.+2)/2 



g£,2A(rf+4)/2(l _Q)3/2 
4*27' 



RE / a;A 



2*2 



ixt2)^fJ,) 



Es 



ADf2 



(r - {xt2Y2^f 

h+t2 [ 2Z?(ii+<2)/3 



xE^p.^ 



N 1 - x9 ^'y EEs , . 2 . N 1 . ^ 

(r - {xt2)^nf - -J—-\(^xt2ll^ - [xi2)^n.v) 



(ti+i2)2'' '"^ (ii+i2) 

(A - l)xt27Vi^+'£^"gKA"r' - A-r') 

47:>3(ti + i2)2A('^+6)/2(l _ 72)3/2(1 _ Q)3/2 



RE ^ xii 

2*7 ilD^ 



-1 - 



Es 



( (r-(xt2)^A)' 

(tl+i2) \ 2D{h+t2)f3 



- 1 



(A - 1)/!'^+^)/^ 
4AL>(l-72) 



,(r-(a;t2)^A)' 



^J^EEs 



_^^2^2_^2(A-1)/} 



(''+')/^ / jEEs ^ ' 

\tl+t2 



{xt2p,'^ - {xt2)^il.v) 



AD 

p(rf+2)/2/_ ... 2 



(/.V-A.r^) 



R= 1 



(A - l)S/|''+"^/"(a;A + t2)E - 2t2Es] + 2Es{xt2)^ (i.v) 



{h+t2){ADA-{A-l)fl 



(rf+2)/2 



xfi^E^) 



(D5) 



(D6) 



(D7) 



(D8) 
(D9) 



A = 1 + 2w/d, and ^^'^ binomial coefRcients. 

Although the above expressions are rather complicated, we are only interested in determining whether the large-d 
behavior in each case is bounded. We notice that all the factors (^2/^1)'^^'^ cancel, and the only terms which have 
a d-dependent exponent are /f^', /2^'i /a^^ and A*^/^. However, the large-d limit of each of these expressions is 
independent of d. We have: 



/f/' = exp(-|) 



hmd-+oo /s^' = exp(- 



2n 

t2X 
2(tl+t2) 



Hmd^oo /2 



d/2 



exp(-f), 
lim^^oo A"''/^ = exp(-w). 



(DIO) 



At large d therefore, all the terms T„ in equation (pO|), are dominated by the exponential factor exp[—dg{x)], and 
hence we may complete the integral using the method of steepest descents. 



2. The 1/d Expansion 

In this section we evaluate the leading order term in the 1/d expansion of each of the expressions for T„ (given by 



equations (|D2|)-(p7D) evaluated at A' = 2D and ip — 7r/2, where ip is the angle between r and A- 



Since the final part of this calculation requires the integration of these terms over the variables x and w, it is 
important to check that these integrals do not alter the c?-dependence of higher orders in the expansion; this ensures 
that we have calculated the entire leading order contribution. If we examine the expression for each r„ in turn, we 
see that, in the large-d limit, every order in the 1/d expansion will have an exponential factor with a negative x and w 
exponent. The presence of these exponential factors ensures that the x and w integrals will not alter the d-dependence 
at any order in the expansion, so we only need to calculate the leading order terms. 
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We first evaluate the variables E, E^, R and Q at the position of the minimum which controls the value of the 
integral in equation (80) {fl — \/2D, tp — n/2); then expanding to leading order in l/d gives E = 1, Eg = 1, R = 1 
and 



w exp(— a;/2) 
^(1-7^) 



X 1- 



(Dll) 



Usi ng t hese results we can now calculate the terms up to 0(1) in the expansion of each expression for r„ (equations 



(D2)-(D7)); these are given by: 



7(^2 — ti) exp{—x/2) w exp(— 1«) 



X U 1 - 



T2 = 
T5 = 



7(^2 ~ ^1)2; exp(— a;/2) w exp(— ti) 
7(^2 ~ ti)x(x — 2) exp(— a::/2) exp(— it)) 
—jr^x exp(— a;/2) exp{—w) 



8i?2(t^+t2)2(l_ ^2)1/2 

-7exp(— x/2) exp(— w) 



815*2 (1 - 7^)3/2 



a; 1 - 



4ti*27^ 

{ti + t2y 



2*27 I-' 



2„2 



D{h + t2y 



t2-h 



t2V^ 



tl+t2 D{ti+t2f 



iti + t2y 



n = o{i/d). 



(D12) 
(D13) 
(D14) 
(D15) 



(D16) 
(D17) 
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